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OF  IHSTABILIUI  OF  STEAM  CCIVECVIVE  IOVEH  EKT  CAUSED  BT  IBTERHAL 
BEAT  SCUBCES 


A.  A.  Takiaov 


The  stability  of  steady  convective  apveaect  caused  by  internal 
beat  sources  was  studied  earlier  in  [1,  2].  Tie  spectra  of  the 
decreaents  of  the  disturbances  and  tie  central  curves  were  obtained 
fpr  different  values  of  the  Prandtl  auabes.  This  report  considers  the 
fen  of  the  disturbances  cf  convective  sovcacnt. 


1.  Heat  sources  with  vcluaetric  dens  it]  Q ere  uniforaly 
distributed  in  a plane  vertical  layer  cf  viscous  fluid  with  a width 
cf  2h d the  layer  is  assuaed  to  be  closed  fro  a above  and  below,  and 
the  vertical  Ralls  which  :bouad  it  are  .held  at  the  sane  tea peratures. 
Stationary  plane- pa rail el  ccnvective  acveaent  with  velocity  and 


i 

[ I 
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temperature  profiles  which  are  even  relative  tc  tbe  axis  of  the 
charnel,  found  frpa  the  ordinary  equations  cf  ccnvection  with 
consideration  of  internal  heat  sources  [1],  originate  in  this  channel 
due  tp  internal  heating.  If  ue  take  fc,  A*/v.  qh*/2,  pg0h*/2 

(q  * a/pcp/)  as  the  units  of  distance,  tine,  velocity,  teaperature, 
and  pressnre,  respectively,  the  diiepsicnless  velocity  and 
tenpexature  profiles  will  le 

v0- 1/60(1 -6xl+5jf4),  (1) 

r0*=  i —x*.  (2)' 

Snail  plane  eprnnl  disturbances 

$(x,  z,  I)  =<r(x)exp(— + (3) 

T (x,  *,*0  -*fl(x)exp(— U + Utz),  (4) 

% 

are  iaposed  on  the  nain  flov,  nhere  ^ and  6 are  the  anplitudes  of 
the  oscillations,  X is  tbe  decrencrt,  acd  k is  the  vave  nunber. 


Substituting  (3)  and  (4)  in  tbe  cesvecticn  equation  and 
considering  the  annllncss  of  tbe  disturbances,  ue  will  obtain  a 
systen  cf  linear  fcoaogeaecus  differential  equations  for  deteraining 
intensities  f (x)  and  d(x)  [ 1 ]; 

A’«f  - ikGHy  + •'  «=  - lA<r,  (5) 

P-'A9  + ikCi  ( Tn  y - vo0)  = - XO  (6) 

(A?  ==*  f"  — kfy,  Hy  ~ vnAy  — v'ny,  G — tfqhr,l 2v*,  P = v/y) 

kith  tbe  hoaogeneons  boundary  conditions 
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(p-«p'-0,  0-0  St  *-±l.  (7) 

2.  le  will  as*  th*  Galerkin  netkod  to  sclvc  boundary  problea 
(5 )- (V)  • He  will  find  cf  and  $ in  th*  fern  (f  the  superinposition  of 
the  tdsic  functions 

X M 

9 = ]£  a, 9b  ,,=!1  biflk-  (®) 

i i *-i  4 

He  will  take  the  intensity  of  th*  disturtacces  in  a quiescent 
fluid*  deterained  fron  the  bcutdaiy  (rctlea 

&9 1 = — \xA9i<  9t  = 9,  = 0,  *t  x = + 1,  (9) 

(/=  1,  2 AT), 

as  the  basic  functions  • end  the  intensity  cf  the  tenperature 
perturbations*  deternined  by  the  pisblct 

P-'AO*  + v*®*  = 0,  \=0  Hx=±l  (10) 

(*=1.2 M) 

as  the  basic  functions  0*  4th*  explicit  fpt«  cf  the  basic  functions 
is  given  in  [3]*  for  •sssfls). 

The  re  gui  resent  of  th*  orthogcnality  cf  the  discrepancies  in  th* 
basic  functions  leads  t*  a syst*a  cf  bciogeteots  linear  algebraic 
equations  for  coef  ficients  s,  and  b„ . 
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The  ccaditioq  of  fche  existence  of  a ncnzccc  solution  for  this 
systen  deter  lines  tbe  spectrua  of  t|e  characteristic  decrements  of 
disturbances  X depending  ca  the  Grashof  ouster  N,  the  Prandtl  nnaber 
f and  the  wave  nuaber  k.  Tic  characteristic  cccieaents  X are  defined 
as  the  intrinsic  values  of  the  systea  aatrix. 


u/  nnn  u ^ 

The  expaqsio*  coef f iciest s^ar 4 the  coa|cnects  of  the  latent 
vector  which  corresponds  to  characteristic  nuafccr  X.  This  vector  was 
icund  as  fellows.  The  valae  of  X was  substituted  in  the 
characteristic  eqaatioq 

\A—KE  I — 0.  (“) 

fche re  A is  the  systen  natriz  end  E is  the  ucit  aatrix.  One  equation 
was  deleted  fron  the  systen  obtained.  Ter  tbe  best  conditionality  of 
the  aatrix  obtained,  tbe  equation  with  the  ainiaua  aodulus  of  the 
coefficient  iq  the  diagonal  tera  was  seleated  as  this  equation.  One 
cf  the  unknowns  (a,  er  b„)  was  assigned  a randca  value  (e.g.,  "*1**), 
aad  the  systea  thes  obtained  (with  ccaplex  cleacnts)  was  solved  by 
the  aethod  of  prinarj  cleaects. 


3.  The  noraal  distorbacces  \f/  (x,  z,  t)  and  T(x,  z,  t)  with  a 
certain  aaplitude  e,  which  teneias  arbitrary  wben  staying  within  the 
beards  cf  the  linear  theory  of  stability,*  ere  ceded  to  the  aain  flow 
aad  distort  it.  fe  will  (let  t|e  current  lines  and  isotheras  of  the 
disturbed  total  aovenent. 


if 
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It  an  arbitrary  fixed  point  it  tiae  tD « the  equation  of  the 
current  line  of  diatarbed  acveaeat  la 

$o(*) +<!*(*,  2,  t0)  *C|, 

khere  0 (x)  is  the  carraat  faaction  of  the  primary  flow,  and  c,  is  a 
certain  constant. 

(f  disturbances  ^(x#  i0  t)  are  deteraired  froa  foraula  (3),  the 
current  line  equation  asaaacs  the  following  fen: 

-V0(.'c)  + al?(x)«,*,  = C,.  (12) 

(actor  eu°  only  a f facta  the  intensity  cf  the  disturbances; 
therefore,  we  can  set  t ( * .0  without  affecting  continuity. 

Considering  t|)e  complex  fera  cf  y (x)  and  ct,  equation  (12)  can 
le  rewritten  as: 

%(x)  + a,  | (x)  cos  -?,(.*)  sin  *z]  = Ax,  (13) 

where 

?,  (■*)  = Re  <f  (x),  f,  (jc)  -=  Im  <f  ( jc), 

ard  a,  and  At  are  xeal  constants. 

Che  isolines  are  plotted  as  fcllonr.  1 be  corresponding  values  of 
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large  ember  of  parts,  as  Mill  obtain  tie  csrrcnt  line  for  the 
selected  value  of  0 1.  Varying  At  with  a certain  spacing,  we  will  find 
the  faaily  of  equidistant  current  lines.  Tie  equation  for  the 
isctheras 

T0(x)  + a,  (0r  (x)cos  Az  — 8,(jr)sln  kz\  = Bx,  (14) 

is  fonad  analogously,  chore 

0,  (x)  = /?«  8 (x),  0,  (x)  = lm  8 (jt), 

and  tie  faaily  of  isotherns  of  the  disturbed  acvenent  is  plotted. 

S.  The  results  of  the  cuaerical  calculaticus  are  given  below. 

The  intrinsic  values  of  tic  systea  aatrix  here  found  using  the  QR 
algcritha  realized  on  Aragats  and  R-22QH  ccaputers  of  the  Coaputer 
Center  cf  the  Perat  state  University  [4],  Six  tc  fifteen  functions  of 
each  type  were  taken  in  expansions  (I) 1 The  latent  vector,  isolines 
and  isotherns  were  found  os  as  it- 2 2 OR  ccnputer. 

figure  1 shows  the  neutral  curve  for  tic  Prandtl  nuaber  P = 10 


plotted  froa  the  nateiials  in  [2].  As  this  study  establishes,  at 
sufficiently  large  values  >ei  t|e  Piardtl  enter,  the  neutral  curve 
ccqsiSts  of  two  branches,  ille  short-wave  branch  corresponds  to 


by dicdynaaic  disturbances  drifting  slcwly  alcng  the  channel.  The 
long-wave  tranch  corresponds  tc  heat  wave  disturbances,  the  phase 
velocity  of  which  is  clcse  to  the  aaxiaua  flew  velocity1. 


fcctnfte:  »Pcr  comparison,  the  broket  liqe  in  Pig.  1 s bout  the 
central  carve  for  t * 0 taken  froa  [1],  which  characterizes  the 
develgpaent  of  hydvodynaaic  disturbances  alcnc.  End  footnote 

It  is  interesting  .to  trace  the  fora  of  the  disturbances  corresponding 
tc  both  branches. 


Eigure  2 shows  the  current  lines  and  isotteras  of  the  total 
coqvective  aoveaent  plqttcd  for  point  A (see  fig.  1),  which  lies  on 
tie  hjdrcdynanic  branch  of  the  neutral  curve  |k  = 3;  G = 2180). 


Figure  3 corresponds  to  pcirt  E,  which  is  located  on  the  thernal 


tranch  (k  * 0.8;  6 = 230).  for  convenience  cf  illustration,  the 
vertical  scale  in  this  figure  is  cne^fgurjbb  cf  the  full  scale. 
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the  values  of  the  curs cat  function  indicated  in  the  fignres  ace 
increased  10*  tines,  and  ti^e  tenperatuce  values  - ten  tines,  then 
cellaring  the  figures  ere  should  reaeaber  that  constant  at  used  in 
toraulae  (13)  and  (14)  a as  2.5  tines  larger  fer  point  A than  for 
point  E then  plotting  the  iselines. 

These  figures  shov  that  in  both  cases,  instability  develops  in 
the  fgra  of  tuo  vortex  chains  ehich  alternate  cn  the  interfaces  of 
the  cpnvective  floys.  Thus,  althcuch  hydrcdjraaic  disturbances  and 
rising  heat  flux  disturbances  are  related  tc  tve  different  nodes  of 
the  instability  spectrua,  there  is  no  cssectial  difference  in  their 
fees.  Hoveveci  the  difference  is  that  heat  saves  have  a relatively 
high  phase  velocity  ccnparcc  tc  bydrodycaaic  disturbances. 
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